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Abstract — Some ceramic materials for high temperature applications are partially transparent for radiative 
transfer. The refractive indices of these materials can be substantially greater than one which influences 
internal radiative emission and reflections. Heat transfer behavior of single and laminated layers has been 
obtained in the literature by numerical solutions of the radiative transfer equations coupled with heat 
conduction and healing at the boundaries by convection and radiation. Two-flux and diffusion methods 
are investigated here to obtain approximate solutions using a simpler formulation than required for exact 
numerical solutions. Isotropic scattering is included. The two-flux method for a single layer yields excellent 
results for gray and two band spectral calculations. The diffusion method yields a good approximation for 
spectral behavior in laminated multiple layers if the overall optical thickness is larger than about ten. A 
hybrid spectral model is developed using the two-flux method in the optically thin hands, and radiative 
diffusion in bands that are optically thick. 


INTRODUCTION 

Some of the ceramic materials being developed and 
evaluated for high temperature applications are par- 
tially transparent to radiant energy. For high tem- 
perature conditions such as in a combustion chamber, 
infrared ant visible radiation from the surroundings 
penetrates into the material heating it internally. Since 
temperatures in the material arc elevated, internal 
radiant emission becomes significant ; this is especially 
true for materials with high refractive indices since 
internal emission depends on the refractive index 
squared. In addition to emission, energy transfer 
within the material depends on internal radiant 
absorption and scattering, and on heat conduction. 
As a result, radiant energy can affect internal tem- 
peratures of some ceramic engine parts and coatings 
that partially transmit radiation in portions of the 
wavelength spectrum. It must be determined when 
radiative processes can be important, and how large 
an effect they have when compared with calculations 
when materials are assumed opaque. 

The refractive indices of single or composite 
materials can have a considerable effect. Surface 
reflections depend on the ratio of refractive indices 
across an interface ; this affects the amount of external 
radiation transmitted into a material, and the amount 
reflected from internal interfaces in a composite. Since 
emission in a material depends on its refractive index 
squared, internal emission can be many times that 
from a blackbody radiating into a vacuum. So that 
radiation leaving through an interface docs not exceed 
that of a blackbody, there is substantial energy reflec- 
tion at the internal surface of the interface, mostly by 


total internal reflection. Scattering is another means 
for energy transfer in the layer; it interacts with the 
internally reflected energy and alters the temperature 
distribution. 

There is an extensive literature on radiative transfer 
in plane layers. Much of the work has been for gases 
with refractive indices of one although some of the 
early work for predicting heat treating and cooling of 
glass plates [1] included refractive index effects. The 
literature has been briefly reviewed in our previous 
work [2-4] and is not repeated here. In refs. [2-4] 
temperature distributions and heat flows in partially 
transmitting materials are predicted by radiative 
analyses using the transfer equations coupled with 
heat conduction. The governing integral equations are 
solved numerically. Each exterior boundary is heated 
by radiation and convection, and the effect of diffuse 
interface reflections is included. Results were obtained 
for a single layer that is either gray or has a two- or 
three-band spectral variation of the absorption 
coefficient, and for a gray Iwo-layer composite. This 
simulates a ceramic layer with reinforcement by 
another material, and also examines the behavior of 
a ceramic coating used to protect another ceramic. 
Various amounts of isotropic scattering are included 
to simulate internal reflections by a granular or rein- 
forcing structure. 

The formulation and solution of the exact spectral 
equations of radiative transfer including scattering is 
rather complicated ; hence it is desirable to have more 
convenient approximate methods if these will yield 
accurate results. Two-flux, diffusion, and hybrid 
methods are investigated here to provide a simplified 
formulation for obtaining approximate solutions. The 
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NOMENCLATURE 


a absorption coefficient of material in layer 
[m '] 

c« velocity of radiant propagation in 
vacuum [ms -1 ] 

C i constant in blackbody spectral flux 

[Wnrsr-*] 

C\ constant in blackbody spectral flux 

[mK] 

integration constant in energy equation 
[W nr 1 ]; CON — CON/DoT 4 g] 

D thicknesses of plane layer [m] 
e vh blackbody spectral flux [W m 2 ] 

FS blackbody fraction in band at small 
frequencies 

G flux quantity defined in equation ( 1 a) 

[Wm '];(?= Gjo T\, 

GS , GL values of G in bands at small and 
large frequencies 

H dimensionless convection-radiation 

parameter, hfaT^ 

h |, h 2 convective heat transfer coefficients at 
two boundaries of layer [W m “ 2 K ' '] 
k thermal conductivity of layer 

[Wm 'K '] 

K extinction coefficient of layer, o + g s [m* 1 ] 
N conduction-radiation parameter, 

k/ffTlD 

n refractive index of layer 

q heat flux [W m : ] ; q - q!<iT\ y 
q+,q radiative heat fluxes in positive and 
negative x directions [W m 2 ] 
q r radiative heat flux in layer [W in -2 ] ; 

q, = qja 7 "*, 

<y r ° externally incident radiation flux 

[W m : ] ; q° — ql'/a 

qS , qL radiative flux in spectral bands with 
small and large frequencies 
R{n) function of refractive index defined in 
equation (22) 

T absolute temperature [K] 


/ dimensionless temperature, TjT^ 

T gl , gas temperatures on two sides of 
layer [K] 

.V coordinate in layer [m] ; X = xjD 
Greek symbols 

k optical coordinate of layer K v 
k 0 optical thickness, K D 

v frequency of radiation [s '] 

p 1 , p° reflectivity of interface for internally or 
externally incident radiation 
g Slefan-Boltzmann constant [W m : K 4 ] 
a s scattering coefficient of layer [m '] 

ft scattering albedo of layer, 

gJ(u + g s ) = aJK. 

Subscripts 
b blackbody 

c value at cutoff frequency 

D based on the length D 

f, s first and second internal interfaces of a 
layer. Fig. 2 

g gas on either side of laminated layer 

h, s higher and smaller refractive indices 

j index indicalingylh layer in a composite 

m the mih frequency band 

r radiative quantity 

S, L spectral bands with small and large 

frequencies 

tot total heat flux by conduction and 

radiation 

v frequency dependent quantity 

1,7 the first and final layers in a composite of 
J layers 

1, 2 quantities at sides x = 0 and D. 
Superscripts 

i, o quantity incident from inside or outside 
-K — flux in positive or negative v direction. 


two-flux equations arc given in refs. [5, 6], and the 
method is shown in refs. [7, 8] to give accurate results 
for gray plane layers with a refractive index of one. 
The method is extended here to include refractive 
indices larger than one, and for heating conditions 
such that the boundary temperatures arc not specified 
and arc found during the solution. The two-flux 
method is found to be convenient for layers of mod- 
erate optical thickness; it yields excellent results for 
gray and two band spectral calculations in single lay- 
ers with or without isotropic scattering. 

The diffusion method is extended here to spectral 
calculations in a multilayered composite where each 
layer can have a different refractive index. It yields a 
very good approximation for spectral heat transfer 


characteristics in single and multiple layers if the over- 
all optical thickness is larger than about ten. 

A hybrid spectral model is developed that uses the 
two-flux method in optically thin spectral bands and 
the diffusion method in bands that are optically thick. 
This also provides good agreement with exact numeri- 
cal solutions. When an optically thick band is present, 
convergence of the iterative solution method is 
improved compared with the two-flux method. 

ANALYSIS 

The two-flux approximation for .spectral radiative 
transfer 

A plane layer of absorbing, emitting, and iso- 
tropically scattering material is subjected on both 
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Fig. I . Geometry, coordinate system, and nomenclature for 
radiative fluxes in the two-flux model for heat transfer in an 
absorbing, emitting, and scattering semitransparent layer. 


boundaries to convection and incident radiation as 
shown in Fig. I . The two-flux method is used to obtain 
the layer temperature distribution and the heat flux 
transferred through the layer. The results are com- 
pared with exact numerical solutions of the radiative 
transfer equations to determine the accuracy of the 
two-flux method for a layer with a refractive index 
greater than one. The boundaries are assumed diffuse. 
The two flux equations are obtained from the radiative 
transfer equations as shown in refs. [5, 6 pp. 772, 
785], and are given by either the Milnc-Eddington or 
the Schusler-Schwarzschild approximations. These give 
the same functional relations, but differ by a numerical 
coefficient in the radiative flux relation. The Milne- 
Eddington approximation is used here. 

The two-flux equations. The two fluxes are in the 
positive and negative directions as shown in Fig. 1 ; 
they are each assumed isotropic and the equations 
are written spectrally in terms of frequency. A flux 
quantity G\ and the net spectral radiative flux q VT are 
related to the two directional fluxes by the spectral 
relations, 

<j\ = 2 (</+ +<y, r ): </„ = </.; ~<h r- (la-b) 

Equations (la) and (lb) can be solved for the positive 
and negative spectral fluxes in terms of G, and </„: 
these relations will be used in the analytical devel- 
opment. 



The two-flux equations including scattering arc given 
in ref. [5] as, 

' d ‘ /r dv = (l-a)14c,h(.v)dv-C r (.v)dv] (.1) 

A, d.v 


1 dC# v 
K, d.v 


dv 


- 3</, r (.v) dv. 


The blackbody spectral flux in equation (3) contains 
the square of the refractive index and is given by 


t\ b (v,r)dv 


27r« : C,v' 
: <T(e rT 


• r -l) 


dv. 


(5) 


A third relation, in addition to equations (3) and 
(4), is the energy equation. For the present conditions 
of steady state without internal heat sources, the total 
heat flow through the layer is a constant. It is given 
by the sum of conduction and radiation heat flows so 
the energy equation is [6, p. 695], 

*--*s + L*- dv - (6) 


The three equations (3), (4) and (6), subject to proper 
boundary conditions, are to be solved for q T (x), G{x) 
and T(x). 

Equation (4) is integrated over all v and the integral 
of q vx over v is eliminated by using equation (6). The 
resulting equation is integrated over .v to yield 

V = -*7-(.v)- £ G -^-dv+CON (7) 


where CON is an integration constant. Evaluating 
equation (7) at .v = 0 and D relates the boundary 
values of T and G to the values of q lM and CON 


CON =kT( 0)+ - 


I pG^I 

U« K 

1 f x G v (7 

q w D= -kT(D)~ ^ 


( 0 ) 


dv 


2^4— dv+ CON. 


( 8 ) 

(9) 


Equation (3) integrated with respect to .v will be used 
later in the form. 


<y vr (x)dv = q n {0)dv+K v (\ -O v ) 


x 


l 


[4c vh (.v) dv — Cr,(.v) dv] d.v 


(10) 


where e vb is a function of x as it depends on the local 
T{x). 

Boundary conditions. The boundary conditions are 
now developed. The total energy flow within the layer 
is by radiation and conduction. Because radiant 
absorption in the two-flux method is a volume 
process, there is no absorption at a boundary surface 
since the surface does not have any volume. Hence at 
each boundary the external convection is balanced 
only by internal heat conduction. The total energy 
flow by radiation and conduction within the layer 
is constant across the layer; the </, ol can then be ex- 
pressed at the boundaries as the sum of external 
convection (which is equal to internal conduction) 
and internal radiation. This yields at x = 0 and D. 


Vu« = A 1 [7; i -7'(0)] + 


J n ¥„(0)d» 


(1 la) 


(4) 
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</,„.= h 2 [nD)-T t2 ] + 


qJD) dv. (11b) l(X) = Jk,! -^.' V '+ r(5 ' V J' < l8 > 


At each boundary the internal radiative flux leaving 
the boundary is equal to the sum of transmitted exter- 
nally incident flux and reflected internally incident 
flux. This yields the relations at x = 0 and Z), 

q: t (0) = q°M-p°) + ‘lM 02a) 

q;AD) = -p°) + q: r (D)p\ (12b) 

Equations (2a, b) are used to eliminate the qf r and 
q~ from equations (12a,b) to obtain expressions for 
G,(0) and G t {D) in terms of the radiative fluxes </, r (0) 
and q yr (D ), 

GAO) = 4 - 2 q„( 0) (13a) 

G X (D) = 4 q° t 2 + 2 z — -- q v AD). (13b) 

\-p 1 -p 

By using equations (8) and (11a), and equations (9) 
and (lib), the T(0) and T(D) are eliminated and the 
following relations are obtained by solving the two 
remaining equations for the integration . constant 
CON, and for ? lot . 



4- q vr ( 0)dv. (15) 

Jo 


Equations in dimensionless form for a gray layer . 
The equations for multiple spectral bands follow quite 
readily from the preceeding relations and from 
relations for a gray layer (one spectral band). The gray 
relations are now given in dimensionless form using 
quantities defined in the Nomenclature. 

Two-flux equations : 

qXX) = ? r (0)+Ml-Q)| [4nV(A')-(?(A')]dA' 

(16) 

G(X) = G(0)-3/c„ | q,(X)dX. (17) 


Equations for G(0) and G(D ) : 


C(0) = 4 1 


( 19a) 

<?(1) = 4 | 


(19b) 


Equations for CON and q Ml : 



Solution procedure by iteration for a gray layer. 
An iterative solution procedure was used by guessing 
values of the dimensionless boundary temperatures 
/(0) and /(l). The boundary relations given by equa- 
tions (19a,b), (18) at X = 0 and 1, and (lla,b) in 
dimensionless form for a gray layer, are solved sim- 
ultaneously to find starling values for q r ( 0) and q r ( I) 
and q lol . A q T ( X) distribution passing through the 
end values and through q {ol at X = 1/2 was used 
as an approximation to start the iteration. Using 
G(0) and G(l) calculated from equations <19a,b) 
the trial q T (X) is normalized such that the condition 
3>K i ^l>q r (X)dX = 5(0)-G(l) from equation (17) is 
satisfied. The G(/0 is then obtained from equation 
( 1 7), and CdN and <y lol are found f rom equations (20) 
and (21). The t{X) is evaluated from equation (18) 
and a new q r {X) is obtained from equation (16). To 
begin a new iteration a damping factor is applied 
between the new and old q r {X) to keep the iterative 
method stable. A small damping factor was needed 
such as 0.0005 when the optical thickness is about 10 
and there is no scattering. Solutions with scattering 
required less damping. Computing time for a solution 
converged to a maximum error on q r {X) of 5 x 10 1 
depends on the integration method used in equations 
(16) and (17). Solutions using a trapezoidal rule took 
seconds on a VAX computer. A Gaussian routine 
look from seconds to minutes depending on the 
optical properties and gave comparable results to the 
trapezoidal rule. 

The surface reflection characteristics were modeled 
by using integrated averages of the Fresnel reflection 
relations. For diffuse incident radiation this gives [6, 

p. US], 


Energy equation : 
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1 (3 n+\)(n-\) 

P (n) =- R(n) = 2 + 6( „ +|) : 


+ 


/r(/r— l) 2 
(n 2 Tf) T 


In 



2n '(n 2 + 2/i — 1 ) 8/i 4 (/i 4 + I) 

(/r+ f)(// J — V) + (/r + l)(/i 4 - 1)- 


x In (/i) (/i = zih/zij. 


( 22 ) 


Equation (22) is for reflection of radiation incident on 
a material of higher refractive index where n h and n s 
are the ‘higher' and ‘smaller’ n values (the incident 
radiation is from within the material with « s ). After 
allowing for energy incident at angles larger than the 
critical angle for total reflection, the p(n) for diffuse 
radiation going in the direction from a higher to a 
smaller refractive index material is found from [9], 


pin) = I — -s[l - R(n)) (n = njnf (23) 
n~ 

Two-flux method for a two-band spectral calculation. 
A two-band calculation is used to illustrate a spectral 
application of the two-flux method. The letters S 
and L designate the ranges with small and large fre- 
quencies. Then for a quantity such as G v (0) the 
band notation is used that <7S(0) = J* c (7,(0) dv and 
GL(0) = J,* G,.(0)dv. For a two-band calculation each 
ofequations(!9) has two parts, one for each frequency 
range. For example, using equation (19a) gives for 
G( 0) in the small and large frequency ranges, 

CS(0) = 4 qS - 2 ff- qSA 0) (24a) 

GL( 0) = 4 ' qL°, - 2 \ + -- : qL,( 0). (24b) 

1 ~P 1 ~P 

Similar relations are written from equation (19b) for 
<5S(\) and GL(l). Equation (21) for q lo{ now contains 
a contribution from each of the two bands. 


quencies. The FS{X) is the fraction of blackbody 
energy in the short frequency range corresponding to 
the temperature t(X)T gl . 

The solution by iteration is similar to that outlined 
for a gray layer except there are now two q t (X ), one for 
each spectral band. The iteration begins by estimating 
<7S r (A') and qL T (X) by taking into account that q r ( X) 
tends to be small in a band where k d is large and that 
the radiative energy distribution in the bands shifts to 
smaller v as t(X) decreases. 


The diffusion approximation for optically thick spectral 
layers 

When the optical thickness of a layer is large the 
radiative diffusion method is very useful for making 
predictions using spectrally dependent properties. 
Since the diffusion method yields equations that are 
convenient to evaluate, the method is extended here 
from a single layer to a multilayer laminated com- 
posite of J layers as shown in Fig. 2. From ref. [6, p. 
753], the radiative flux by diffusion in any of the layers 
is given by d q fT = — (4/3AT v ) (de vh /d.v) dv. To obtain 
the radiative heat flux including all frequencies this is 
integrated over all v to give 


H r = - 




(28) 


For integrating over all v it is convenient to use mul- 
tiple spectral bands where K v = K„, in the wth band 
and there is a total of M bands. Then with each band 
extending from v„, to v„,+ , where 1 ^ m ^ M, 


Hr 



e vb (v)dv. 


(29) 


The quantity in the integral can be expressed in terms 
of the blackbody fraction F(v/c 0 T ), 



*Vb(v)dv = n 2 oT 4 lF(v m + 


Jc { J)-F(vJc,T)l 


Hun = H 


<75(0) GL( 0) 


CON 

N 


+ qS r (0) + qL r {0). (25) 

This also occurs in the temperature relation equation 
(18), 


(30) 


The total energy transfer by radiation and con- 
duction is then given by. 




dT 4 
k d.x + 3. 


I 

K m 


t(X) = 



GS(X) 


GL(X) 

3/v"i>i 


-q m X+C0N 


X dv I // 2 <T T" 4 !^ v„ 


,lc 0 T)-F(vJc u r 1| . (31) 


i 

» • 


(26) 

Equation (16) for the radiative flux is written for each 
band. For the band with short frequencies. 

qS r {X) = qS,(0) + K m (\-n.) 

x f [4/rV(Ar)+S(AT)-GS(A')]d4r (27) 
and similarly for c/L r (X) in the band with large fre- 


Sincc q iol is a constant through the composite layer 
equation (3 1 ) can be integrated from 0 to x / in the yth 
layer. The result is placed in dimensionless form using 
the quantities defined in the Nomenclature to yield, 

H\o\X, = N t [t /f — t f (X,)]+ ^ 

m - 1 
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Fig. 2. Laminated multilayered geometry showing coordinate system and nomenclature designating 
interfaces and their temperatures for use in spectral diffusion method for single or multiple layered 

composites. 


— F{vJc 0 tj f T g] )} — t*(Xj){F[v m + \lc 0 tj(Xj)T g] ] 

-F[vJc 0 t,(X t )T^)) (32) 

where t jS = t, (X t = 0) is at the first interface in the y'th 
layer (see Fig. 2). 

Equation (32) is evaluated at X i = 1 which gives an 
expression for the temperature difference across the 
yth layer in terms of q xov To join two layers a con- 
tinuous temperature is assumed at the internal inter- 
face between them. The temperature jump that occurs 
in some instances at an interface when using the 
diffusion approximation is neglected at the bound- 
aries. The jump is small when the layers are optically 
thick and when conduction is comparable to radiation 
[10]. The convective and radiative balance at each 
outside boundary, and the conduction with radiative 
diffusion relation for q iox , equation (32), written for 
the temperature difference across each of the J layers, 
provide the following simultaneous equations: 

//i(l-/. f ) + 0~P i .>i[^-^]-^o, = 0 (33a) 

4/7 2 

M | 

x L 7 n -['J ,lc 0 i lf T e[ ) 

m - I 

~ F(v,Jc tt t /f T t] )} — t* s {F[v m+ i/c 0 //. 5 7’ g ,] 

- F[vJc 0 t hS T g} ]}] (1<7<7) (33b) 

tf2(0.s-'„2)+(l * 0. 

(33c) 

These equations are solved for the two external 
boundary temperatures / Kf and /,- s , the internal inter- 
face temperatures, and the heat flow q loi through the 
composite layer. After <y lol is obtained, the temperature 
distribution in each layer (1 ^ j ^ J) is calculated by 
solving equation (32) numerically for t,{X f ). 


Hybrid method using the two-flux method for thin 
bands and the diffusion method for thick bands 

The hybrid energy equation. In the numerical solu- 
tions using the two-flux method, convergence was 
found more difficult for a layer with optical thickness 
larger than about 20. To analyze spectral cases with 
both optically thin and thick bands, a hybrid method 
was devised. For optically thin bands the two-flux 
method is used, while the diffusion method is used for 
thick bands. For simplicity a two-band model is given 
here where the optically thick region is for 0 ^ v ^ v c 
so that k os > k dl . The total flux by conduction and 
radiation is 



Using the diffusion relation from equation (28) for 
the 1st integral and the two-flux relation, equation (4), 
for the 2nd integral gives 

d T(x) 4 d f r * 

~k £f c - (r,d> ' < 35 > 

By integrating equation (35) and putting it in dimen- 
sionless form the temperature distribution becomes 

1 f 4 

/(*) = - v— t(X) 4 FS[t(X)] 

N ( 3 k ds 

GL(X ) } 

— — —QwiX+CONf (36) 

as compared with equation (26) using the two flux 
method for both spectral bands. Since t(X) in equation 
(36) is in the blackbody fraction /\S[f(A')], a numerical 
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method such as a root solver is used to find the dis- 
tribution t(X). 

Boundary conditions. For the diffusion approxi- 
mation, external radiative energy in the optically thick 
band is absorbed and emitted at the boundary surface 
at ,r = 0. The net external energy that is absorbed, 
combined with the external convection, is equal at the 
boundary to the internal heat conduction combined 
with the internal radiative diffusion (which acts like 
conduction) in the optically thick band (small v). The 
total internal heat flow at the boundary can then be 
expressed as this net external energy (absorbed radi- 
ation and convection) added to the internal radiation 
in the optically thin band (large v) which is obtained 
from the two-flux method, 

<7u, - (1 -<jT(0) 4 FS[T(0)]} 

+h l lT 9l -T(0)]+ql40). (37a) 

Similarly at .v = Z), 

9,„, = -(\-pl{qS"^aT(DrFS{T(D)]) 

+ h 2 [T(D)-T i2 ] + < i L,(D). (37b) 

From equation (36) in dimensional form evaluated 
at .v = 0 and D, 


CON = kW)+ — ■ 7’(0) 4 FS[7'(0)]+ ~ GL(0) 

jA s 3a l 


(38a) 


Viol 


CON k 
~D D 


T(D) 


- vlD TiD) ‘ FS ' T(D>i - rKj> c ^ < 38b > 


Combine equations (37a) and (38a) to eliminate the 
T(0) term, and combine equations (37b) and (38b) to 
eliminate the T(D) term. Then eliminate q lo{ from the 
resulting two equations. This yields the equation for 
CON in dimensionless form (as in equation (20)) as 


CON = 



'( l+ 


GZ.(1) 4 , 1 

+ + ... t( I ) F5[/( I )] 


3jVky 


(0) + (l-p°) 


x [qS^ - t(0) 4 FS[tm) J —qL r ( I ) + ( 1 - pi 

x [</£?: “ d I ) 4 FS{?{ l )]]J • (39) 

The q xox can be obtained from from either equation 
(37) or (38b). 

Solution method. The iterative solution starts by 
assuming qL{X) s /(O) and t{D). The G£(0) is found 
from equation (24b) and similarly for GL(1). The 


t]L r ( X) is normalized according to the two-flux equa- 
tion (17) so that 3K nL y o qL T (X)dX = GL(0)-GL(\). 
Then GL{X) is found from equation (17) as 

GL(X ) = GL{0)~ 3k dl ^ qL r ( X) dX. (40) 

The CON is obtained from equation (39), and q xox 
from equations (37) or (38b). The t(X) is then obtained 
from equation (36), and a new qL r (X) is computed 
from the two-flux equation (27) written for the large 
v band. To start the next iteration a damping factor 
is applied between this qL r (X ) and the qL r ( X) that was 
used to start the iteration. 


RESULTS AND DISCUSSION 

Tw o- flux results for a gray layer 

The most basic application of the two-flux method 
is for a gray layer. Two-flux results are given in refs. 
[7, 8] for a layer with a refractive index n = I between 
two black walls at specified temperatures. Excellent 
agreement was obtained with numerical solutions in 
the literature of the exact radiative transfer equations. 
In the present work the boundary temperatures are 
unknown and are obtained in the solution as a result 
of specified convection and incident radiation con- 
ditions at the boundaries. The present results show 
the effect of having a refractive index as large as n = 4 
and the effect of isotropic scattering with an albedo 
of ft = 0.9. 

Figure 3 is for n = 2, and Fig. 3(a) is without scat- 
tering. The ‘exact’ numerical results were obtained 
in refs. [2, 3] where the behavior of the temperature 
distributions was discussed. The two-flux method 
yields excellent predictions of the temperature dis- 
tributions and the results follow the pronounced tem- 
perature curvature near the boundaries. There are 
only small deviations from the exact solution; the 
largest deviations are for k* d = 1 where radiation has 
a large effect since the layer is neither optically thin or 
thick. The total heat flux through the layer by com- 
bined radiation and conduction is predicted very well 
by the two-flux method as shown on the figure. 

The effect of large scattering, ft = 0.9, in Fig. 3(b) 
is for the same n and K n values as in Fig. 3(a). For 
the same k 0 a layer with ft = 0.9 has 10% of the 
absorption as for ft = 0. For these conditions the two- 
flux relations provide excellent predictions for all 
of both the temperature distribution and total energy 
transfer through the layer 

When the refractive index is increased to n = 4, with 
all other parameters kept the same as in Figs. 3(a),(b), 
the results are in Figs. 4(a), (b). The exact numerical 
solutions are again obtained from refs. [2, 3]. The tem- 
perature profiles are more uniform in the central por- 
tion of the layer as discussed in refs. [2, 3] which leads 
to large temperature variations near the boundaries. 
The two-flux method is successful in predicting both 
the details of the temperature distributions and total 
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Dimensionless position, X=x/D Dimensionless position, X-x/D 



0 . 3 * 1 1 * 1 1 1 

0 . 0 0 . 2 0 . 4 0 . 6 0 . 8 1.0 


Dimensionless position, X=x/D 

Fig. 3. Comparisons of two-flux results with exact numerical 
solutions of radiative transfer equations without and with 
scattering in a gray layer with refractive index n - 2; < 7/1 
= 1, <ft = 0.25\ - 1 , = 0.25, //,-//.- l, N = 0.1. 

(a) Scattering albedo. 0 = 0. (b) Scattering albedo, O = 0.9. 



0 . 3 l . 1 * ' * » L 1 1 

0 . 0 0 . 2 0 . 4 0 . 6 0 . 8 1.0 

Dimensionless position, X=x/D 


Fig. 4. Comparisons of two-flux results with exact numerical 
solutions of radiative transfer equations without and with 
scattering in a gray layer with refractive index n = 4 ; c/m 
- 1, = 0.25\ / gl = 1. t t2 - 0.25, H x = H, = I, N = 0.1. 

(a) Scattering albedo, H = 0. (b) Scattering albedo, Q = 0.9. 


heat flux transferred by combined radiation and con- 
duction. 

Two-flux results for a layer with two absorption bands 
For the results in Fig. 5(a) the spectral absorption 
coefficient has two values : in the small v range (vJc a T $l 
< 1/4000), a iS D= 10, while in the large v range 
(vJcoT $ \ > 1/4000), a tl D = 0.1. The selected value of 
the cutoff frequency v c divides the blackbody spectrum 
at T g] into approximately two equal parts. A gray 
scattering component is present so that for the two 


temperature profiles in Fig. 5(a), o^D — 0.9 or 
a S D = 9.9. Thus for the two curves the spectral optical 
thicknesses are, respectively, k ds = 10.9. k, >l = 1, and 
k* os = 19.9, k ol = 10. As shown the two-flux method 
applied spectrally, as developed in the analysis, pro- 
vides excellent agreement with the two-band numeri- 
cal solutions in ref. [3], 

In Fig. 5(b) the a v D values arc increased to 100 and 
I for the bands with small and large frequencies, and 
the gray scattering component is o s D - 9. The spectral 
optical thicknesses are then k !>s = 109 and k 1>l = 10. 
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Dimensionless position, X-x/D Dimensionless position, X=x/D 



0. 0 0. 2 0. 4 0. 6 0. 8 1.0 

Dimensionless position, X=x/D 

Fig. 5. Comparisons of two-flux and hybrid results with 
exact numerical solutions of radiative transfer equations for 
two-band spectral calculations with scattering and a refrac- 
tive index of n - 2; q' r \ = 1. q l ,\ = 0.25 J , ( ti = I, = 0.25. 
H x = = I, N = 0 1, vjc lt T^ { = 1/4000. (a) Band absorp- 

tion thicknesses : a^D - 10 <a vl D - 0.1. (b) Band absorption 
thicknesses: a yS D = 100, a vl D = 1. 

The two-flux method provides excellent agreement of 
t(X) and </ u „ with the exact numerical solution of 
the transfer equations from ref. [3]. Results from the 
hybrid method are also shown and they provide very 
good agreement. For a large in one band the hybrid 
method provides more rapid convergence than the 
two-flux method. 

Spectral diffusion results for single or multiple layers 
Temperature distributions are shown in Fig. 6 for 
a layer with refractive index n — 2 and three spectral 
bandi> within which the optical thicknesses are 5, 10, 


Fig. 6. Comparison of spectral diffusion method predictions 
with temperature profiles from exact numerical solution of 
the radiative transfer equations for a layer with three spectral 
bands without scattering and with an albedo of 0.9 ; q° y = 1 , 
q? 2 = 0.25 4 , , = 1 , /,: = 0.25, H x = H z = 1 , jV = 0. 1 . 

Optical thicknesses in the three bands in order of increasing 
frequency are k d , — 5, 10, 20, and the two cutoff frequencies 
are given by v c /r 0 r g , = 1/5300, t/3300. 

and 20 in order of increasing frequency. The two 
cutoff frequencies (given in the caption) divide black- 
body radiation at T %y into approximately three equal 
parts. Two exact numerical solutions are given : the 
solid line is without scattering (Q = 0) while the long 
dash line is for Q — 0.9. The diffusion approximation 
is independent of the amount of scattering and hence 
only one diffusion curve is given (short dashed line). 
The diffusion curve is close to the other curves but 
does not follow their exact shape near the boundaries. 
This is expected as the diffusion approximation (based 
on an isotropic assumption) does not apply adjacent 
to a boundary because the flux in this region is not 
isotropic. The diffusion method predicts the heat 
transfer through the layer within 10%. 

The diffusion method provides a very good pre- 
diction of the temperature distribution in a three layer 
region as shown in Fig. 7. The exact numerical solu- 
tion was obtained by extending the computer program 
from ref. [4]. Each layer has three spectral bands 
defined in the caption. The diffusion results are con- 
siderably easier to compute than the numerical solu- 
tion of the radiative transfer equations for three 
coupled layers and three absorption bands. The 
opaque limit for conduction (no radiation) within the 
layers, and with emission, absorption, and convection 
at the outer boundaries is shown as a dot-dash line. 
This shows that there is a substantial internal radi- 
ation effect although the combined optical thicknesses 
through the three layers is moderately large. A good 
prediction (within 4%) is obtained for the heal trans- 
fer through the composite. 
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Fig. 7. Comparison of spectral diffusion method predictions 
with temperature profiles from the exact numerical solution 
of the radiative transfer equations for a three-layer composite 
with three spectral bands; $?, = 1, q? 2 = 0.25 4 , r gt = I, 
t f2 = 0.25, Hi = H 2 = I. Optical thicknesses in the three 
spectral bands in order of increasing frequency are 
K Di , = 5, 10, 20, K 0y = 10, 5, 5, k D j , = 5, 20, 10, and the two 
cutoff frequencies are given by vJc 0 T gl - 1/8000, 1/4000. 



FiG. 8. Comparison of hybrid two-flux/diffusion method 
with the exact numerical solution for a layer with one opt- 
ically thin spectral band {k dl = 1), and one optically thick 
band (K m =100); ft = 1, ft 2 = 0.25 4 , /,,*!, t t2 = 0.25, 
//, = H 2 - 1, N = 0A. The cutoff frequency is given by 
v c /c 0 r Rl = 1/4000. 


Results using hybrid method for spectral calculation 
with one hand optically thick 
The two-flux solution computed by iteration had 
more rapid convergence when the diffusion method 
was used in the optically thick band rather than the 
two-flux method. Results from the hybrid method for 
a two-band calculation are in Fig. 8 for refractive 
indices of n = 1 and 2. The boundary temperatures 
are predicted quite well and the hybrid curves deviate 
only moderately from the exact solutions within the 
layer in the regions where the curvature is largest. 
Good predictions are obtained for the heat transfer 
q xoi through the layer. 


CONCLUSIONS 

The prediction of temperature distributions and 
heat transfer is carried out in semitransparent layers 
heated on both sides by radiation and convection so 
that the boundary temperatures are not specified and 
must be determined in the solution. The layers have 
refractive indices larger than one, and isotropic scat- 
tering is included. Three methods were developed for 
performing spectral calculations. The two-flux 
method was found to give excellent agreement with 
exact numerical solutions for all the conditions con- 
sidered. The diffusion method was developed for spec- 
tral calculations in multilayered media and was found 
to give good predictions for large optical thicknesses. 
The diffusion results are much easer to evaluate than 
a numerical solution of the radiative transfer equa- 
tions for a multilayered region with several spectral 
bands. A hybrid method was developed that can be 
applied when there are optically thin and optically 
thick spectral bands. This gave good temperature and 
heat flux predictions, and was easier to evaluate than 
the two-flux method when a spectral band is optically 
thick. 
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